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Abstract 

If stable Heavy neutrinos of 4th generation possess their own Coulomb- like interaction, recombination 
of pairs of Heavy neutrinos and antineutrinos can play important role in their cosmological evolution 
and lead to observable consequences. In particular, effect of this new interaction in the annihilation of 
neutrino-antineutrino pairs can account for 7-flux observed by EGRET. 

1 Introduction 

This work begins systematic study of model of subdominant component of dark matter in form of Heavy 
neutrinos pQ in the special case when this component possesses its own long range interaction. It is supposed 
that new interaction is Coulomb-like, being described with unbroken J7(l)-gauge group. We call it y- 
interaction. Its massless gauge boson and charge are called y-photon and y-charge, respectively. The 
existence of this interaction as well as of 4th generation (and 4th neutrinos), possessing it, can naturally 
follow from superstring inspired phenomenology 2 . It is assumed according to £Q that Heavy neutrino 
(N) belongs to new (4th) generation of fermions, that the electroweak and strong charges of 4th generation 
quarks and leptons are attributed analogously to other three Standard Model generations and that mass of 
4th neutrino lies in range about m — 45 — 80 GeV. Strict conservation of y-charge implies the lightest fermion 
of 4th generation (4th neutrino) to be absolutely stable. Due to y-charge 4th neutrino cannot have Majorana 
mass, while its Dirac mass is of the order of other fermions of 4th generation. Thus y-interaction provides 
natural basis for the hypothesis of massive stable 4th neutrino, y-charge neutrality implies neutrinos to be 
accompanied by the same amount of antineutrinos. 

In the early Universe effects of new interaction can influence the freezing out of y-charged 4th neutri- 
nos. The existence of y-photons implies y-radiation background (y-background), interacting with y-charged 
neutrinos. 

After decoupling from plasma and y-background 4th neutrinos and antineutrinos due to their long-range 
y-interaction can form bound systems. These bound systems annihilate rather quickly. Therefore such 
"recombination" of 4th neutrinos and antineutrinos (N-recombination) reduces their relic density and can 
lead to observable effects. 



In the Galaxy, 4th neutrino pairs, being strongly subdominant component of Dark Matter (DM), can 
enter the clumps, formed by the dominant DM component and the enhancement of N-recombination in such 
clumps can lead to 7-fiux, saturating high energy diffuse 7-background, observed by EGRET. 

The existence of sufficiently longliving y-charged quarks of 4th generation would add new element 
to the evolution of 4th neutrinos and the self-consistent picture for cosmological multi-component y-plasma 
requires separate consideration. Therefore we consider below the evolution and observable effects of 4th 
neutrinos and y-background only. 



2 y-plasma and y-radiation in early Universe 

In the early Universe at the temperature 1 T > m 4th neutrino pairs were in equilibrium with cosmological 
plasma, what lead to exponential decrease of their equilibrium concentration at T < to. When the annihila- 
tion rate became less, than the rate of expansion (at Tf ~ to/30) 4th neutrino pairs were frozen out and their 
amount in the co-moving volume did not change essentially in the course of subsequent expansion. How- 
ever, at lower temperatures, after decoupling of the frozen out y-charged 4th neutrinos from y-background, 
formation of bound systems neutrino- antincutrino (N-recombination) was possible. 

Annihilation timescale of systems with the size cib bound due to interaction with the fine-structure 
constant a y is given by 

r ~ ^ (1) 

Or 

V 

and it is quite short. Therefore such recombination is in fact equivalent to effective annihilation of neutrino 
pairs and it should reduce the number density of primordial 4th neutrinos. 

Let us first consider the influence of y- interaction on the process of freezing out (at T ~ m/30) of 
4th neutrinos. New interaction opens new channel of 4th neutrinos annihilation, which is analogous to 
27-annihilation of e + e~ 

NN -> yy. (2) 

Furthermore for slow 4th neutrinos and antineutrinos (with relative velocity v < 2ira y ) it leads to the increase 
of cross section of all their annihilation channels, including one via intermediate Z-boson (NN — > Z — ► ...), 
due to Coulomb-like factor of Sakharov enhancement 

C = rAv C = 2 -^. (3) 

1 — e °° v 

At neutrino mass of interest (to ~ 50 GeV) channel NN — ► Z — ► ... is close to resonance (m res — mz/2 w 
46 GeV). Therefore the effect of channel ©, being suppressed near m res , grows with the deviation of to 
from m res . 

The ratio of cross sections of annihilation channels of (J2J and through Z-boson in non-relativistic limit 
is independent of Coulomb factor and given by 

^ W 144f-gg-V Pz[m) (4) 
az \l/30j P z (50GeV)' 1 ' 



^2 \ 2 r 2_2 

Pz(m) = 1 



Am 2 J 16to 4 ' 

where mz and Tz are the mass and the width of Z-boson, (50 GeV) ~ 0.0289. Effect of Coulomb factor 
also does not lead to the principle change, since at T ~ m/30 velocity of 4th neutrinos is still large. 
Number density of 4th neutrinos n can be expressed through its ratio r to the entropy density s 



r 



45 



T'\ (5) 



where g s takes into account all the species of ultrarelativistic bosons and fermions and T is the temperature 
of photons. 

1 Throughout in the text the system of units H = c = k = 1, where k is the Boltzman constant, is used. 
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Figure 1: Number densities of 4th neutrinos frozen out at T ~ m/30 in units of entropy density for the cases 
without and with y- interaction with different magnitudes of its constant. 



On the figure 1 the frozen out number density of 4th neutrinos (r = — ) is presented. 
Both effects of channel (01 and Sakharov enhancement induce correction of density of 4th neutrinos 
within the factor 2-10. 

After freezing out the gas of 4th neutrino pairs is for some period in thermal equilibrium with ambient 
ordinary matter and y-photon background. 

y-background has direct thermal exchange with 4th neutrino pairs only and due to their small frozen 
out concentration it can not experience essential thermal exchange and be in thermal equilibrium with 
other matter (ordinary plasma and radiation) at T < Tf. In the course of successive evolution the number 
density of y-photons reduces relative to that of photons by a factor k < 1. Due to entropy conservation, the 
less species remain in equilibrium with radiation, the smaller is the value of k. In the period of Big bang 
nucleosynthesis (at T ~ 1 MeV) k rs 1/6.5, and y-radiation can not influence signficantly the predictions for 
light element abundance. 

4th neutrinos interact with each other due to long-range y-interaction. This interaction with a large, 
Rutherford-like cross section provides their equilibrium (Maxwell) distribution, whatever energy losses or/and 
exchanges with other matter components they experience. 

Before decoupling of 4th neutrinos from y-radiation their temperature Tv was equal to that of y-photons, 
i.e. Tjv = k 1 / 3 T. After decoupling of 4th neutrinos from y-background (at T <C Tv y ) their temperature 

T 2 rr.. .oouv( ™ \ 3/2 V30 



T 



Ny 



" V50GeV/ a,. 



The probability for 4th neutrinos and antineutrinos to form their bound states grows when particles slow 
down, so at T < Tv y the rate of this process increases. For estimation of this effect we will use classical 
approximation following 4 . TV and N moving towards each other due to y-attraction must experience dipolc 
emission of y-radiation. If this radiative energy loss exceeds the initial energy of their relative motion, they 
become bound. The typical length for such energy loss is much less than the mean distance between 4th 
neutrinos (a fortiori than Debye radius of N-plasma), and the timescale for this energy loss is much less than 
the timescale of N-N and N-y interactions. It provides such binding to proceed freely, independent of other 
interaction processes. Cross section of bound system formation in this approach [5] is obtained to be 

cr b =irpl(v)= J u/B d , (7) 

where is the maximal impact parameter at which NN pair is bound due to y-radiative energy loss, v is 
the initial NN relative velocity and a c i = a y /m. 

The rate of binding with cross section having the form ok = ctq/v^ for particles, distributed by Maxwell, 
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IS 

4r(2-f) /m\^ 
r rec = ( <J b v )n = 2 p ' a I — I rs. (8) 

After N — y decoupling T rec exceeds the expansion rate H at T below 

In this expression ry ec denotes the value of r(m, a y ) before recombination starts, being determined by fig.l, 
and r(50GeV,l/30) = 1.22 • 1CT 15 . 

The most of created bound systems have typical size less than pf, ~ a c i/v 7 ^ 5 ~ ac^m/TV) 7 / 10 ; and their 
annihilation timescale Eq.Q turns out to be much less than the timescale of their destruction as well as 
than cosmological timescale. 

During recombination, 4th neutrino gas will be effectively heated, since slow pairs recombine more effec- 
tively and disappear. From the other hand, due to dipole emission of NN pairs scattered at p > p b without 
binding, NN-ga.s cools down. Evolution of 4th neutrino pairs with the account for the both thermal effects 
can be described by the system of equations 

dT = \ a b v /HT 

Here variable 6 shows deviation of neutrino temperature from Eq.© due to thermal effects, E and E re i are 
the kinetic energy of a single neutrino (or antineutrino) and of relative motion in NN pair, respectively. 

3 1 3 

^2 Tn ~ E ~ 3^ I ' ei ) f7bV ) = ifl^bV^N, (11) 

where 7 = 5/ \"^ n = \- Note that the terms ((§TV - E)a b v) and ( {-\E re i)a b v ) in Eas. l|10llljl describe 
thermal effects due to pair disappearance (heating at (3 > 1) and due to dipole emission (cooling), respectively, 
and the first one prevails (7 > 0). The system (|1U|) can be solved by iterations (dividing the first equation 
by the second equation one gets r(9) = r rec # -1 / 7 and then subsitutes this value into the system again). It 
gives for r at RD stage (at T > T R m ~ 1 eV) 

r= — r- (12) 

{l + F R (T,T rec )} lh 



and at MD stage (at T < T RM ) 



{1 + F R (T RM , T rec ) + F M (T, T RM )Y h ' 

The functions F R m have the form 



(13) 



F Rt M(T,T ) = r recJ ^-D RM - -^j . (14) 



Here a = - 2 = §, M = - § = ±, 7 = 1 + = §, 

D R = DmV5W = (a b v) {dT2) ^ S « 1.6 • 10 18 eVt f * ( SSl)* . (15) 

v H K 1 H R T V50GeV/ \1/S0J v ; 

Hubble constants on RD and MD stages are related as 



where g £ takes into account contribution into energy density of all ultrarelativistic species. At T <C To the 
term l/Tg R > M in Eq.iO can be neglected. Also 1 + F R can be neglected in Eq. l|13[l . Note that without the 
account for thermal effects (0=1) 7 = 1. 
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Figure 2: Relic densities of 4th neutrinos in units of critical density for the cases with recombination, without 
it and without y- interaction at all. a y = 1/30, 1/60 were taken. 



In the period of galaxy formation gas of 4th neutrino pairs becomes strongly nonhomogeneous. Since the 
contribution of 4th neutrinos into the total density is negligible, they cannot play dynamically important role 
and follow the dominant component of dark matter in the process of development of gravitational instability. 
That is why evolution of 4th neutrinos in this period strongly depends on the model of galaxy formation. 
According to CDM model of galaxy formation at z ~ 10—20 the distribution of Dark matter becomes 
strongly non-homogeneous at small scales. Such structure can not appear in the model of Hot Dark Matter, 
which is however strongly disfavoured by the observational data. 

In the course of development of gravitational instability the fraction of particles, which remain outside the 
structure of nonhomogeneities decreases with time. Therefore N-recombination in homogeneous gas of 4th 
neutrino turns out to be suppressed, when galaxies began to form at T = 7/j„ = 10T mo( j, where T mo( i = 2.7 
K. Recombination of 4th neutrinos in Galaxy with the account for small scale nonhomogeneity (dumpiness) 
of their distribution is considered in Section 4. 

On Figure 2 the relic densities of 4th neutrinos are shown in units of critical density for a y = 1 /30 and 
1/60. For comparison such densities without effects of recombination and without y-interaction at all are 
put on Fig. 2. 

If (meta-)stable quarks of 4th generation (Q) exist and possess y-charge, they represent another compo- 
nent of y-plasma, and joint evolution of y-charged N and Q should be studied. However, preliminary analysis 
indicates that there is a range of Q-gas parameters, at which presence of Q-component does not influence 
significantly the main features of N-recombination in early Universe. 



3 7-emission from recombination of primordial 4th neutrinos 

Among the products of N-recombination (and annihilation) on RD stage only ordinary (light) neutrinos can 
survive to the present time. Their maximal modern energy, being red-shifted from the period of such early 
N-recombination, does not exceed few MeV. It makes their relatively small flux hardly detectable in neutrino 
observatories. 

Hadronic and electromagnetic cascades induced by N-recombination interact with the ordinary plasma 
and radiation. The energy, realized in such interaction leads to distortion of CMB spectrum. If energy 
release in a unit volume exceeded 10 _4 e 7 and took place later than at T w 5 keV [H], it would be observed 
in CMB spectral measurements. Total energy realized in the result of recombination at T < 5 keV is given 

by 

Se 2ms(5keV)r(5keV) 8 (h keV\ * / r rec \ A /50 GeV\ * /1/30 

e^ < £ 7 (5keV) ~ \~T~) V1.22- 10- 15 / \ m ) 

This estimation makes annihilation effects be hardly constrained by the data on CMB. 
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Figure 3: Gamma fluxes from recombination of 4th neutrinos in early Universe in comparison with EGRET 
data. a y — 1/30, 1/60, 1/137 and m = 50, 75 GeV were taken. 



If N-recombination proceeds on MD stage, 7-radiation induced by annihilation in this period will con- 
tribute into extragalactic 7-emission in the energy range measured by EGRET. 

During N-recombination, within interval, when Universe temperature falls down on dT, s ■ dr pairs of 
NN in unit volume annihilate, where dr is given by Ea. (|10|) . 

Let Bz = 77 — ^ — be fraction of annihilation acts going through the channel NN — > Z — > being 
determined with the help of Eq.JSJ), and dN 1 (E ) — f-y(E )dEQ be averaged multiplicity of created 7 in this 
channel within energy interval E Q . . . E + dE Q . Due to redshift modern energy E . . . E + dE of photons 
emitted at the temperature T (redshift z) is given by 

E = (z + l)E=-^-E, (17) 

^ mod 

Taking into account that number density evolves as oc s, for modern number density of 7 products of early 
NN annihilation we have 

dn^E) = B z ■ dN^Eo) ■ s(T mod ) ■ dr(T) = B z s mo d ■ f-y(E {T, E))-^—dE ■ r' T dT, (18) 

mod 

where r' T = s mo d — s{T moc i) fa 3 • 10 3 cm~ 3 . One can pass from integration dT to dE . For intensity 
from Ea. (|18|) we have 

^ {E) = = BzS Z d E2 m ° d Cj" A(30 • MT{E, Eo)) ■ E dE . (19) 

T mo d E 

Here Eo max is the upper limit of annihilation 7-spectrum, c is the light speed. The relationship between Eq 7 
E, T is given by Ea.lfTT)). 

Using results of previous section one obtains intensity of 7-radiation from primordial 4th neutrino re- 
combination given on Figure 3 for different values of m and a y . EGRET data on extragalactic 7-radiation 
are shown for comparison. Spectra f 1 {Eo) had been obtained with the help of code PYTHIA 6.2. 

Effect of TV variation (growth) due to disappearance of slow pairs and to dipole emission was taken into 
account here. This effect raises the predicted 7-flux by factor 2-3. 

The growth of Tm reduces ( er&u ) , what seem to contradict the mentioned increase in the predicted 7-flux. 
The fact is that we are interested in annihilation effects produced on the late stages of N-recombination. 
Indeed, since the maximal energy of the created 7 is equal to the mass of 4th neutrino, annihilation photons, 
born in period T < 2000 K < T rec only, contribute to the energy range of EGRET. To this period the bulk 
of the primordial 4th neutrinos had been annihilated, so the net effect in the interested range is proportional 
integrally not to the initial r rec but to their residual number density, which is in the inverse dependence on 
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( <JbV ) . By the same reason the inverse dependence of predicted 7-flux on a y is explained, being explicitly 

22 3 

given by I 7 (a y ) oc r' T (a y )B z (a y ) oc a y 23 r?£ c (a y )B z (a y ). 

Comparison with EGRET data restricts the allowed values of a y as a y > 1/137 for m — 50 GeV. This 
constraint becomes slightly more severe for larger m. 

Note that the products of N-recombination at T < 2000 K can provide early ionisation of neutral matter 
gas, which is favoured by WMAP data. 



4 Effects of y- charged 4th neutrinos in Galaxy 

Due to their negligible contribution into the cosmological density (see Fig. 2), 4th neutrino cannot play any 
significant role in the dynamics of galaxy formation. Their distribution in Galaxy will be governed by the 
dominant DM component, which we assume to be CDM, for definiteness. In the gravitational potential 
of Galaxy 4th neutrinos acquire large virial velocities u ~ 200 km/s. It strongly reduces the rate of their 
recombination so that the corresponding timescale highly exceeds the age of the Universe t mo d- However, 
the typical lifetime of created bound systems in Galaxy is much less than t mo d and slow N-recombination 
can produce 7-flux, accessible to observations. 

Assume that 4th neutrinos are distributed in Galaxy according to isothermal halo model with density 
profile 

(n\ (lk P c) 2 + (8.5kpc) 2 

p(R)=PNloc (lkpc)2 + fl2 , (20) 

where R is the distance to the galactic center, pmioc = £ ■ 0.3 GeV/cm 3 . The distribution lj2U|l is determined 
by the dominant CDM component and the 4th neutrino density is taken to follow it, being proportional to 
the ratio £ of 4th neutrino and dominant CDM densities. Velocity distribution is assumed to be Maxwcllian 
with the maximal likelihood value u = 220 km/s. It turns out that the 7-flux, predicted for various halo 
models is rather weakly sensitive to their choice, varying within a factor 1-2. 

The account for small scale inhomogeneity (clumps), of dominating Cold Dark Matter (CDM) strongly 
enhances the predicted 7-flux from N-recombination. The possibility of such a small scale structure and its 
survival to the present time is determined by primordial perturbation spectrum and on dynamics of galaxy 
formation |7| . We will take into account by parameter rj the enhancement of N-recombination (annihilation) 
rate due to dumpiness. In the case of 4th neutrinos such enhancement can be much larger, than for non- 
interacting DM particles, since 4th neutrino annihilation in clumps goes through N-recombination, which 
is enhanced not only due to increase of concentration, but also (and dominantly) owing to much smaller 
virial velocities of 4th neutrinos inside the clump, as compared with their averaged density and velocity 
distributions in the Galaxy. 

The following effects are important in the estimation of parameter 77. 

Since clumps are formed by the dominating CDM particles, the minimal mass of such clumps is deter- 
mined by the ultraviolet cut in the spectrum of their density fluctuations, which is determined for weakly 
interacting particles by the scale of their free streaming. As in pQ we'll use the results of [7] in the de- 
scription of such CDM clumps. The existence of the period of thermal equilibrium with ordinary matter 
and y-radiation makes gas of 4th neutrino hotter, than ordinary CDM gas. It prevents 4th neutrinos from 
entering the smallest CDM clumps and there exists the minimal mass of CDM clump, in which the presence 
of 4th neutrinos is not suppressed. This critical minimal mass is estimated to be equal to 

„ 2» r /50GeV\ 15/4 / a v \ 3/2 

where Mq is the solar mass. Following we can suppose that 4th neutrino density distribution at scales 
larger than M^min is proportional to that of CDM, so their relative contribution is £ = £In/&Cdm, where 
i7jv is given by Fig. 2, Qcdm — 0.3. Radius of clump with M = M/v mm is equal to 
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Typical velocity of particles inside the clump with M = Mn m i n is determined by 



GM N 

u c i = d ^ wl9m/s. (23) 

We will assume that this velocity parameter is constant over all the volume of clump. Here we do not consider 
the possibility for 4th neutrinos to enter the clumps of mass less than Mm m i n , which was considered in 1 . 

Clumps of mass around the minimal one give the main contribution into the annihilation enhancement 
effect. The shape of the number density distribution n inside the clump determines the effect of annihilation 
by the factor 

v r , 



S =JpJ n ^ (24) 

where V — ^nR 3 , N = J ndV is the total amount of given particles inside the clump; n = ^ gives the mean 
number density. 

In case of 4th neutrinos, annihilation rate strongly increases due to its velocity dependence so their 
number in clump can experience essential change. In our study of evolution of 4th neutrinos in the clump 
we will take into account effects of heating due to slow pairs disappearance and of dissipation due to dipole 
emission (the latter turns to be less important, than the former). Parameter Eg. (|23l) is connected with the 

2 

initial temperature of N: T/vo = — ^ ■ Any decrease of the total number of 4th neutrinos in the clump 
does not affect its gravitational potential and mass, being determined by dominant DM component. We will 
assume that the shape of distribution of 4th neutrinos inside the clump does not alter with time, and that 
the value of the factor l|24l) for 4th neutrinos can be taken the same as for non-interacting DM, which was 
estimated in |7| to be S « 5. Then for total amount of 4th neutrinos inside the clump one has 



dN _ I _ „. \ a t-2 S 

— _- {<JbV )N v ^ 



-{(T b v)N 2 

l^f = { (l T N -E- \E rel )a b v )iV|. 



Solution of this system is analogous to that of Eas. (|10[) and yields 



' (26) 



Here N ann means the total annihilation rate given by first equation of system l|25|l . Index "0" corresponds 
to the initial moment T/v = Tjvo- The mean number density of 4th neutrinos in the clump and characteristic 
time of recombination in it r are estimated as 



no 



1.8 • 10" 8 cm- 3 



50 GcV 



2.310- 



0.20 GyrM^ (_»*_) I (3g>) * 



(27) 



where £(m = 50 GeV, a y = 1/30) = 2.3 • 10~ 8 was used; note that £(50,1/137) = 1.1 • 10~ 7 . Neglecting the 
dependence of r on r rec , one obtains £ cx m 18 / 23 a y 22 / 23 ; so T i s almost independent of m and a y . 

The obtained value of r exceeds the clump dynamical timescale ~ R/ud w 0.005 Gyr what provides 
the shape of 4th neutrinos distribution to follow the "equilibrium" form inside the clump. However, the 
temperature Tjy during t — to^ 10 Gyr becomes 2 times larger. It implies that for the "equilibrium" form of 
N distribution the effect of N evaporation can be important. It can lead to the slowing down of annihilation 
rate in the clump as compared with the above estimation. 

Enhancement factor r\ for 4th neutrinos, corresponding to the present moment may be estimated as 



9 / 5 / ^ \ *+i , _ . q/4 / -i /on\ 9/10 



{—J < 3 - 10 *l»G5v) 8/4 (^) ' (28) 

Here the factor t]q oc nci/ n Gai 0- e - proportional to the ratio of averaged squared number densities of DM 
particles inside the clumps and steeply distributed in Galaxy) denotes the enhancement of annihilation of 
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Figure 4: 7-fluxes from recombination of 4th neutrinos in early Universe and in Galaxy with enhancement 
factor 200 in comparison with EGRET data for a y = 1/137 and to — 50 GeV. 



non-interacting DM in clumps. From figures 3-5 of [7] one obtains its value in the range 1 < 7/0 < 5 for 

Mjnin — Mn min • 

The estimation (128(1 does not take into account the mass spectrum of clumps. The latter is predicted 
to start from the minimal value M = M m i n and to fall down with the increase of M. The account for this 
mass distribution of clumps will slightly decrease the value of 77. Indeed, in Eg. (|28|l the u c i dependence on 
M gives additional suppression ex M~ 3//5 . 

Active annihilation of 4th neutrinos inside the clumps after they separate from the cosmological expansion 
at z ~ 10 — 20 contributes into the intergalactic 7-background. However, in the period of clump formation, 
when the density inside the clumps is of the same order as the average density (Sp/p ~ 1), this contri- 
bution cannot exceed the effect of N-recombination by homogeneously distributed 4th neutrinos. During 
the successive periods of galaxy formation the amount of neutrinos inside the clumps decreases due to their 
annihilation, while homogeneously distributed neutrinos enter the inhomogeneities and their fraction reduces 
with time. It leads to the suppression of extra-galactic contribution into the high energy part of gamma 
background as compared with the local effect of N-recombination in clumps of the halo of our Galaxy. 

Intensity of the modern 7 radiation from annihilation of 4th neutrinos in the direction I to halo of Galaxy 
is determined by 

m = V .^flf l{ E)[{£-) 2 d l. (29) 

For to — 50 GeV, a y — 1/137 with factor 7/ w 200 one obtains 7-flux shown on Figure 4. 

To saturate the high energy part of EGRET data at larger values of a y , a slightly larger factor rj is 
required. For instance, if a y = 1/30 (at m = 50 GeV), 77 « 240. Such a weak sensitivity of the predicted 
galactic 7-fl.ux to the value of a y follows from the fact that this flux depends on a y virtually only through 
such dependence of Bz ■ The latter is weak for the values of to close Z-boson resonance. It strongly differs 
from the case of primordial recombination, where in addition neutrino velocity depends on a y through T/v y . 
That is why a slight change of the factor 77 can easily compensate the decrease of Bz{a y ) for larger a y . 

5 Discussion 

The existence of new long range interaction would lead to new forms of matter, such as y-plasma and to a 
set of new phenomena, in which N-recombination, considered in the present paper plays important role. 

In the considered case of y-charged 4th neutrinos N-recombination suppresses strongly the pre-galactic 
density of this subdominant component of Dark Matter, but appears possible to explain all the data on 
extragalactic 7-emission observed by EGRET by the effect of N-recombination for to 50 GeV and a y w 
1/100. 
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It should be noted that in the framework of Grand Unification (GUT) models the value of a y can hardly 
be less than a em — 1/137. Moreover, such a small value of a y implies rather specific requirements for GUT 
model, which embeds y-interaction. Indeed, y-interaction, being possessed by 4th generation fermions only, 
should be unified with other interactions at large energy scale. Therefore one can expect that it should have 
constant at low energies greater than electromagnetic one, since the last is reduced from the unified (GUT) 
constant due to screening by virtual particles of all the four generations. 

The future theoretical model, which will embed the considered here new long range U(l) interaction, 
should satisfy a set of necessary conditions, such as anomaly freedom of y-interaction. The y-charge assign- 
ment e v of 4th generation particles (N, E, U, D) like e V N = e y E = —e y u/3 = —£yD/3 would avoid triangle 
anomalies within Standard Model field content extended to four generations. However, triangle anomaly 
problem, as well as the possible effects of y-7, y-Z mixings |Sj should be considered within a unified frame- 
work. To avoid theoretical inconsistencies such framework can involve possible extension of particle content 
(SUSY partners, exotics of Eq superstring inspired GUT model, etc.), what requires separate study. 

Taking into account these theoretical problems, the present analysis has no aim to give any final conclu- 
sion. It just opens the room for the systematic search for WIMPs, possessing new long range forces. One 
can apply to these searches the methods, developed in this work. One can hardly imagine the impact of the 
existence of such particles. Reminding Zeldovich, we can only repeat his words: "Though the probability 
for these phenomena is small, their expectation value is great!" 
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